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Abstract 

We investigate the signatures, in the squeezing limit of the primordial scalar bispectrum, of 
modifications of the standard theory at high energy, such as a modified quantum vacuum or modified 
dispersion relations, reflecting the effects of different physics at higher energies/earlier times. Our 
results are obtained through the analysis in field theory and not using the (approximated) templates 
for the bispectrum. In this way we can obtain the precise dependence on the probed scale and 
quantify the presence of enhancements and/or non-local shape of the non-Gaussianitics, which 
are relevant, for example, for measurements of the halo bias. It is shown that the signatures are 
distinctive and can be determined in general, and the results are different from those obtained using 
the templates proposed in the literature, where existing, in particular because of the interplay 
between the scale of the observations and the scale of the high-energy physics when taking the 
large-scale (squeezed) limit. Several pieces of information regarding high energy physics could be 
obtained in case of detection of these signals, especially bounds on the scales of new physics. 
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A Appendices 

1 Introduction 

The non-Gaussian features of the primordial cosmological perturbations are going to become one 
of the most important experimental and theoretical tools for investigating the early history of the 
Universe pp. The scalar bispectrum, coming from the three-point correlator of the comoving curvature 
perturbation £, is the leading contribution to the scalar non-Gaussianities in the perturbative theory. 
In recent times, several reasons of interests have emerged for the investigation of its so-called squeezed 
limit, when one of the three momenta it depends on in Fourier space is much smaller than the other 
two: kx < ks, k 2> 3 ^ k s . 

In fact, large-scale probes of the primordial non-Gaussianities like the halo bias are expected to 
be soon competing with or even surpass in accuracy those of the CMBR and this prompts for a 
better understanding of their sensitivity on the details of the bispectrum and the theoretical models. 

Furthermore, studying the squeezed limit provides powerful results capable of distinguishing (hope- 
fully falsify) classes of models. For example, all single-field models (adiabatic perturbations) with 
Bunch-Davies vacuum and standard kinetic terms (we call this the standard scenario from now on) 
predict the so-called local form of the bispectrum and a very low level of non-Gaussianities in the 
limit [3HB]. Remarkably, measurements of the halo bias could discriminate cases with different under- 
lying assumptions 0[8]: for example models with modified kinetic terms whose bispectrum template 
is maximized on equilateral configurations, or with initial state different than the Bunch-Davies, in 
which case the proposed template is maximized on folded configurations [9] . 

In all the analysis in the literature, however, the halo bias and the squeezing limit have been 
investigated using the templates and not the actual field theoretical results for the bispectrum 0. 
Sometimes the discrepancies between templates and field theoretical results disappear in the squeezing 
limit: for example, Maldacena's result [I] becomes factorizable in that limit, and assumes indeed a 

1 While completing this paper, we have found that in [TO] the analysis and the simulations regarding the halo bias are 
still done using the templates, but the authors warn that they could lead to the wrong predictions for the observables. 
However, the detailed study of the field theoretical results was not done, in particular not for the folded template. After 
the appearance of this article, we were also signalled QT|, which also makes some remarks, focusing on the presence 
of enhancements, about the squeezing limit in a particular formalism (density matrix approach) with a specific initial 
condition for the field mode functions. However, also in this case the detailed study of the general field theoretical results 
was not undertaken (the conclusions in that paper are in agreement with those in one of the classes of scenarios we deal 
with). 
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local shape. However, we will show that in other cases the use of templates in place of the actual 
theoretical predictions appreciably changes the outcomes. 

In this paper, we want to look at the field theory cosmological predictions for the squeezed limit 
of the bispectrum for single-field slow-roll inflationary models, investigating the effects of two very 
simple modifications of the standard scenario: 1) the presence of higher derivative corrections to the 
action that are unsuppressed at a certain physical scale A and modify the dispersion relations of the 
perturbation fields [12], and 2) the choice of an initial state other than the Bunch-Davies vacuum, 
reflecting the effects of different physics at higher energies/earlier times as proposed in |13y i4j. 

The description of the physics at high energies is largely unknown, and, at the same time, there 
are various different models of single-field inflation. It is therefore important to derive general re- 
sults, as indeed was in the spirit of We will therefore deal with this topic adopting the same 
phenomenological approach of [9l ll5t fT6]. not assuming specific detailed high-energy models, but still 
obtaining quantitative as well as qualitative results. We will present in more details the approach and 
the scenarios in section (3[ It is also important to stress that we consider the realistic squeezing limit, 
where k\ is small but non-zero, as defined in [6HH], in view of possible future observations. 

A part from [9lll5j . examples of the modification to the vacuum that we describe have been proposed 
in |13U14| and many successive papers. On the other hand, modified dispersion relations occur in a 
variety of models, such as certain realization of quantum gravity, as the Hofava-Lifshitz model [17j . 
or also in braneworld models [18], and are of interest for the studies of Lorentz violation [19j . They 
occur as well in the effective theory when the fields propagate with small speed of sound |20j. 

Before presenting the detailed analysis, it is useful to outline the points that could make the 
difference with the standard scenario in the cases we consider. In fact, the standard result follows from 
a simple argument [3HS]: in the standard scenario, the leading contribution to the non-Gaussianities 
comes from late times, and furthermore, in the squeezing limit, we consider one of the three momenta 
(k\) that the bispectrum depends upon to be very small. Therefore the leading contribution to non- 
Gaussianities occurs when k\rj <C 1 (r] is the conformal time), and thus the perturbation depending 
on k\ is at superhorizon scales, is frozen and acts as a background for the other perturbations in the 
bispectrum. It has the effect of shifting the horizon-crossing time for the other perturbations, and 
the result is thus determined by the (slow-roll suppressed) scalar tilt of the spectrum, see [IHS]. It 
has been argued that, because of the generality of this argument, this result for the squeezing limit of 
the bispectrum is valid for any single-field slow-roll model -with Bunch-Davies vacuum and standard 
dispersion relations-, with all possible couplings that allow a slow-roll inflation. 

However, this argument would fail for models with an additional important contribution to the 
non-Gaussianities at much earlier times when it still is k\r\ > 1 and the perturbation is not frozen at 
superhorizon scales. The scenarios of modified initial vacuum and modified dispersion relations have 
in fact such additional early time contribution to the non-Gaussianities: it is due to particle creation, 
respectively at the time when the condition picking up the non-Bunch-Davies vacuum is imposed, 
and at the time/scale of appearance of the new physics that modifies the dispersion relations. One 
could think that the particle creation would represent a tiny effect, being severely constrained by 
backreaction (as we will review), but the bispectrum can be quite sensitive to it, because particle 
creation will lead also to interference effect and phase cancellations in the correlators that could lead 
to enhancements. 

We want to investigate the effects of these contributions on the squeezing limit of the bispectrum 
by performing the analysis at the rigorous level of the field theory description, investigating whether 
general features appear and if violations of the standard result occur. It is important to make two 
comments regarding our work. First, the analysis we make is general, within the scenarios of section 
[3l and ranges over all possible single-field slow roll models, allowing all those operators that permit a 
slow-roll inflationary evolution. Although many operators would be quite constrained or suppressed 
because of backreaction/preservation of slow-roll, as it was discussed in [6[l21j . the allowed ones are still 
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numerous, and one could think that looking for the emergence of general features would be hopeless. 

However, the squeezing limit allows powerful results that rely on very general features of the 
scenarios, common to all the specific models, permitting, as we said before, to constrain entire classes 
of models. In the case of the standard scenario, the constraints following from the freeze-out argument 
we reviewed above are so strong that the full leading form of the squeezing limit of the bispectrum 
is obtained, valid for any model of that scenario, while the differences among the specific models 
appear only in the subleading (by powers of ^) corrections. Also in the modified scenarios we are 
considering, there are general constraining arguments (particle creation, interference phenomena), and 
we will show that the squeezing limit is in fact dominated by these general aspects of the scenarios, 
at leading order, and thus there appear general features in the squeezed bispectrum. 

Second, let us recall that in the case of modified initial vacuum, one can divide the contributions 
to non-Gaussianities in two parts [9l[T5]: one deriving from intrinsic non-Gaussianities at the bound- 
ary where the boundary condition defining the initial state is imposed, the other from the particle 
production due to the modified vacuum. The first kind of contribution is absent for Gaussian initial 
state modifications. As in [9j[13j[15], we focus on this class of modifications, which represents the 
minimal deviation from the Bunch-Davies vacuum. We will show that already this minimal change of 
initial conditions has quite an effect on the squeezing limit. Moreover, at least in the boundary action 
formalism, the leading intrinsic non-Gaussian modifications have been shown to be of the local type 
-see |22j and the comment in [9]- and, being also quite constrained by backreaction, they would thus 
lead to a squeezing limit similar to that in the standard scenario. 

In the case of modifications of the initial quantum state, we will also show that the field theoretic 
result in the squeezing limit is very different from those obtained in previous studies, see 0[S], using 
the folded template proposed for CMBR analysis [9]. This disagreement could have been already 
anticipated, since the standard evaluators (cosine and fudge factor) for the matching between the 
proposed template [9j and the theoretical prediction [15] generically indicate that the two depart more 
and more for large kLT] c , where ki is the largest momentum and r\ c is the time when the boundary 
condition picking up the non-Bunch-Davies vacuum is imposed, see [9]. Even more importantly, the 
template does not depend on the parameter rj c , and thus taking the squeezing limit in the full result 
might (and will) be different that doing it in the template, because of the presence of distinct scales. 

The outline of the paper is as follows: we start in section [2] with a very rapid presentation of the 
necessary formalism and notation, followed by a brief review in section of the scenarios that we deal 
with. The analysis of the realistic squeezed limit of the bispectrum is reported in sections 0] and [5j 
First, we present in sections 14.11 and 14.21 the detailed computations for two examples of interactions 
(which will be also proven to yield among the largest contributions to the bispectrum), to illustrate 
step-by-step how their contributions are affected by the features of the modified scenarios. Then, we 
generalize the study to the other various possible couplings in section 14.31 The analysis and discussion 
on the results are contained in section [5j Finally, we briefly comment on the consequences of our 
results on the halo bias in section [6l and conclude in \7\ Some useful equations and material, which 
would have made heavy the main text, have been moved to the appendices. 



2 Formalism and notation. 

We begin by introducing the basic notation and elements of the formalism of cosmological perturbation 
theory, for more details see [2]. Scalar perturbations in single-field slow-roll inflationary models are 
efficiently parametrized by the comoving curvature perturbation^ 

—3(k(vy k - x . (i) 

(27r) 2 



2 We follow the conventions in [H ll5l[To1 ; other authors call this perturbation TZ to distinguish it from the uniform 
energy-density curvature perturbation. 
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The slow-roll parameters are defined as e = „ 2 f.^ , 77 , = —-^7 + » „ 2 jfx > H is the Hubble 

rate, (j) the inflaton (background), and a(rj) ~ — jj- is the scale factor at zeroth order in slow-roll. 
Throughout the paper, dots (primes) indicate derivatives with respect to cosmic time t (conformal 
time 77). 

The two-point function is defined as 

(C(kiK(k 2 )) = (27rf5^(k 1 + k 2 )P(k), k = \h\ . (2) 
and, following [23], we define the bispectrum B as 

<Ck>)Ckfa)Cfc 8 fa)> = (27rf5(^h)B(k 1 ,k 2 ,k 3 , V ). (3) 

i 

A key ingredient in the computation of correlators is the Whightman function 

(2tt)W& + k 2 ) Gfarf) = (CfefoXfeO/)) = i^f^Hh + k 2 ) (4) 
where f^(v)jft(v) are two linearly independent solutions of the field equation 



/£+H^) 2 -y)4 = ° 



2 z \ „ a 

such that, once quantized 



C^) = ^\ + JL ^H- (6) 

Imposing the standard commutation relation on the operators cL ^cL, entails a certain normalization 
for the Wronskian of /|(?7)j /^( 7 ?)- Different choices for fkiv) are different choices of vacua [3]. 

In equation {5]), u is the (comoving) frequency as read from the effective action: we will limit ourselves 
to the isotropic case, where uj depends only on k, dropping the arrow symbol. 

3 Brief review of the scenarios under consideration. 

To ease the discussion of the results for the three-point function and make the paper self-contained, 
let us recap here the main features of the scenarios we will deal with, for more details see [91 I12I - H6] . 

3.1 The standard scenario 

In the standard scenario the dispersion relation is the standard Lorentzian one and the mode functions 
ftiv) are determined by choosing the Bunch-Davies vacuum, so that 

k) 2 = k 2 , f k (v) = ^e^+ l f ^H^i-kv) , (7) 

where is the first Hankel's functions with v = § + 1 ~ 2 " s , and the equations and solutions are 
extrapolated up to 7/; n = —00 where the initial state has been chosen as the empty vacuum (Bunch- 
Davies). n s = 1 — 6e + 2t7 s i is the spectral index. For small momenta, at leading order in e, 

Planck 

where the suffix st indicates that we are considering the standard scenario. 
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3.2 Modified dispersion relations 

Higher-derivative corrections to the action are deeply rooted in the well-established effective field 
theory framework. In fact, one expects that at certain times at least some of the higher derivative 
corrections to the action cannot be treated perturbatively, because, since the metric scale factor a{rf) 
is rapidly decreasing back in time, the suppressing factors, given by powers of ratios such as ^ = a A A , 
where A is the relevant high-energy scale, are not small any more. These no longer small corrections 
can lead to modified dispersion relations. As recalled in the introduction, such possibility appears in 
a variety of scenarios, see for example those in [17H20j . 

As we said, we adopt a phenomenological approach and consider generic modifications to the 
dispersion relation due to higher-derivative corrections to the effective action, encoding them in a 
function F{— ■?■) by writing the comoving frequency as 

Ufa k) = a(ri)Lo phys (p) = a( V )p F (-£) = k F (jkrA , F(x ^ 0) ^ 1 , f » 1 • (9) 

In [16}l24j. it has been shown that if the modified dispersion relations always satisfy the WKB 
conditions (adiabaticity) at early times (that is, WKB is violated only when the dispersion is in the 
standard linear regime u ~ k as usual) the bispectrum is similar to that of the standard scenario, and 
so it will be its squeezing limit. This is easily understood, since the particle production is very limited 
in those cases. 

Instead, [16] has shown that the bispectrum is particularly sensitive to modifications that lead to 
a violation of the WKB conditions (adiabaticity) at early times, where particle production is more 
substantial. We will therefore focus on those scenarios. 

The generic shape of a dispersion relation with WKB violation once at early times is depicted in 
figure [TJ We stress that we are not proposing any specific model, which would be a strong assumption 
on the high energy physics: we consider the most generic function F (— x) i n equation Q, with the 
only constraint that adiabaticity is once violated at early times (it is easy to extend our analysis to 
multiple violations). We are in fact interested in the general consequences on the bispectrum from this 
scenario. Similar phenomenological approach and generic dispersion relations have been considered, 
studying the spectrum, for example in the references [12], see also [19j . 

A general treatment of the field equation (0) in these cases yields the following solution |16| : 



fk(v) = { 



ftui(77, k) I : 7] < rjr 

B l Uxirj, k) + B 2 U 2 (ri, k) II : 4 k) <V< vi k) 



a T* n i( 7 ?' k ) + PT* u 2(?/, k) III : rj[ k) <r] < r)^ 

Dx Vifa, k) + D 2 V 2 (v, k) IV : r$ < r, 
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where ui,2(?7> k) can be very well approximated using the WKB method, while Vi.2(?7, k) and U\p.{?)-, k) 
need to be found by other means, as the WKB conditions are not satisfied!!. As known, given equation 
([5]), the times i/ I I ]n where the WKB approximation first fails are in proximity of the turning points 
of V(rj, k) 2 = lo(j], k) 2 — — , see |25|I26|. They depend on k; however, when clear from the context, we 
will omit the label ( k \ The coefficients D\ : 2, Bi^,o^ dT , (3^ dT are obtained by asking for the continuity 
of the function and its first derivative, and by imposing the Wronskian condition W{/, /*} = — i in 
order to have the standard commutation relations in the quantum theory. 

The details of the solution can be found in [16], and are briefly reviewed in the appendix I A. 11 
Some of the salient features needed here are that: 

• backreaction constraints the interval of WKB violation to be very small, that is 

Vi 

expanding for small A, we obtain generically (i.e., independently of the detailed form of u;(r], k)) 



mdr / .V^fyrjr Q 
a k = 1-* 7T 



V 2 



VI 



V 2 



(12) 

where V(rj, k) 2 = u)(ri,k) 2 — Here, y%\ i is the order 1 factor signalling the violation of the 
WKB conditions at 77,, and Q is reported in equation (jlOOp in appendix IA.ll The magnitude 
of |/9 fe r | is constrained by backreaction. We will discuss more of the ^-dependence and the 
magnitude of these coefficients in section 13.41 

we choose ft = 1 picking up the usual adiabatic vacuum. 



mdr 



The two-point function in the small momentum limit, at leading order in e, (5 k , A, is 

P QA El ( l + 2Re(/?r dr A • (13) 

Planck v ' 



3.3 Modified initial state 

Let us review the two implementations of the modified vacuum (initial state) approach that have been 
first proposed in |13yi4| , reflecting the effects of different physics at higher energies /earlier times than 
inflation in setting the initial/boundary condition for the perturbation fields: 

• EBA approach: fix the boundary conditions for the fields at the same specific finite time r\ c (be- 
ginning of inflation) independently of the modes k, through an effective boundary action (EBA) 
which accounts for the high energy physics via renormalization |14j . This kind of boundary 
condition strongly breaks scale invariance, being imposed for all modes k at the same bound- 
ary time. The (cutoff) scale A, proper of the effective action formulation, is not related to the 
boundary time r/ c . 

• NPHS approach: consider the effective theory valid up to a certain energy scale A, and choose 
the boundary conditions for the solution to the field equation at the new physics hyper-surface 
(NPHS) corresponding to when the physical momentum reaches that cutoff scale |13j . In this 
case, the cutoff is imposed in a scale invariant way via the condition ; = A, so that the time 

rj c when the initial state is picked is k- and A- dependent. At zero slow-roll order, r\ c = —-j^jj- 

3 Some different general approximations are possible, see [IB] and the appendix lA.il With particular choices of time 
coordinate, also the WKB approximation can be applied sometimes, see |25| . 
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To spare notation, we use the symbol rj c for both the NPHS and EBA cases and make clear in the 
context which approach we are following. 

As in [HldS], when studying these scenarios we further adopt a phenomenological method and 
consider generic Gaussian boundary conditions. The latter is a minimal change from the Bunch- 
Davies vacuum. In this case the three- and higher- point functions are determined by the two-point 
function. However, recall that in principle the boundary condition defining the modified vacuum can be 
chosen to be intrinsically non-Gaussian, leading to further contributions to the three-point function not 
determined by the two-point one. Let us also point out that, in fact, in the EBA language, it was shown 
in [22J -see also the comment in [5]-, that the leading corrections due to intrinsic non-Gaussianities in 
the boundary condition would be of the local form, which would thus yield the same form of squeezing 
limit as in the standard scenario (with a suppressed amplitude, because of backreaction and lack of 
cumulation with time, see later sections). We will however focus on the minimal Gaussian vacuum 
modification, as done in [SlUS] (see also the study of the related observational prediction's in [7HS]) 

Then, both in the EBA and NPHS cases, the solution of the field equation is of the Bogoljubov 
form: 

f k ( V ) = arV^vHW(-k V ) + PFV^HPH-kri) , (14) 

where a™ v , /3™ v are determined by the specific boundary conditions (choice of vacuum) imposed on the 
solution at the time r] c , and by the Wronskian condition, which translates into (a™ v ) 2 — (/3™ v ) 2 = 1- in 
general the Bogoljubov coefficients depend on the mode k and the cutoff scale and time (these latter 
are related in the NPHS case, but not in the EBA, as we said). 

The two-point function at small scales, at leading order in slow-roll and , is 

Planck v / 

3.4 General constraints on Bogoljubov coefficients 

The precise magnitude and scale (k) dependence of the coefficients /3™ dr ,/3™ v , entering respectively 
equations (|13p for the modified dispersion relations scenario, and (|15p for the modified vacuum, are 
determined in full details by the specific model or boundary condition giving rise to the modified 
dispersion relation or modified vacuum. However, to be in accordance with observations and to avoid 
backreaction stopping the slow-roll inflationary evolution, these coefficients are subject to a series of 
phenomenological constraints, see for example [15|I16|: 

• the observations on the spectral index show that — fc f^2i£jL^l*i)) i s small, so that /3™ dr , /3™ v must 
be slowly varying with k at the observed scales; 

• backreaction imposes two further constraints: 

— the total energy density must be finite. In general this demands that |/3™ dr |, |/?™ v | decay 
faster than k~ 2 at large k; 

— preserving the slow-roll inflationary evolution enforces the constraint |15|.I16J 

W%Wn < y/M HM £°* , /* = UU - e • (16) 

In the following we will keep indicating the scale dependence of these coefficients with the label 
observe also that we distinguish the modified vacuum scenario from that with modified dispersion 
relations with the labels "mv" and "mdr" . 
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4 Bispectrum in the squeezing limit 



We will now discuss the three-point function. We start by presenting two examples of calculations 
for two representative interactions (with and without higher-derivatives), to show in details and step- 
by-step how their contributions to the bispectrum are affected by the generic features we outlined in 
the introduction (particle creation, presence of new scales related to new physics and/or initial time). 
These interactions were studied, but not in the squeezing limit, in [Hl9|[T5 | [T6 | l2T] • 

We will then generalize our analysis to the other possible cubic interactions and see if the above 
mentioned generic features of the scenarios suffice to determine in general, and in what measure, the 
leading order form of the bispectrum. Finally, we comment on the results in section 

Let us introduce some notation. For the purpose of computing the bispectrum, it is convenient to 
use the field redefinition jll9l [T5fl 

C-C + - ^ C 2 + - ^ d~ 2 (C d 2 C ) + -^C 2 (17) 

Planck Planck T 11 

Using the definitions of the slow-roll parameters, the bispectrum of £, which is the relevant one for 
observations, is related to the three-point function for £ c as [I] 



2k? 



» i=l 
mod 3 

The two-point function and the quadratic part of the action are the same for £ and Q c [3]. 

Using the variable £ c , the general formula for the three-point function at leading order (tree-level) 
in perturbation theory is 

(Cc(v, XiXc(?7, %2)(c(v, £3)) = -2Re (^J drj'i(iJji n \C, c (r], Xi)Cc(»7, x 2 )( c (r], £3 )#(/)(?/) I , (19) 

where r/i n , \ipm) are the initial time and state (vacuum), and Fm is the interaction Hamiltonian. 

The modifications due to the scenarios we are considering will affect both the first connected 
and time- integrated contribution on the right-hand side of equation (|18p . both the other not time- 
integrated ones, via the modified two-point functions. However, the latter modifications, as we have 
reviewed in the previous sections, are strongly suppressed. We will thus focus on the time-integrated 
connected contribution, which is expected to be more sensitive to the modifications appearing in our 
scenarios, see [T5j[16]. 

We will write the result for the three-point function in powers of |/3fcJ. Up to linear order, 

(C, k, C c> g 2 C Cj fc 3 > = So (C Ci ^ C c> k 2 C Cj fc 3 ) + Sp (C C) ^ C Cj fc 2 C c , k 3 ) (20) 

with 

So(C c fc 2 C k 3 ) = A (h,h, k 3 )5F (k u k 2 , k 3 ) , 5 P {C C ^C c g 2 C c k 3 ) = A(h,k 2 , k 3 )SF 1 (k 1 , k 2 , k 3 ) 

(21) 

In the standard scenario we define <Ji^ tandard = i ; ^standard _ so that A(ki,k2,ks) is indeed the 
standard three-point function for Q c . In the modified scenarios, 5Fq and 5F\ represent therefore the 
corrections to the standard result, respectively of order 1/3^1° and due to the new high energy 

physics. 



4 In [15] . ^"c is simply written as £, see their equation (3.11). The same is done in many other papers, such as [9]. 
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4.1 Example 1: Minimal coupling cubic interaction 

We begin by studying the case of the cubic interaction [5] : 

H {1) = ~ [ c*W - JL_ g8r*£ . (22) 

J Planck 

We call this interaction the "minimal coupling cubic interaction" because it is already present in the 
simplest case of a scalar field (inflaton) minimally coupled to gravity [3]. 
In this case equation (|19p leads to 

<C £l (V)C C , k 2 (V)C C> n 3 fa)> =2Re ^(2tt) 3 S(Y^ h) (j^j f dr/^- f[ d v ,G ki (r?, 7/)+permutations) , 

(23) 

where r/ ~ is a very late time when all modes ki are outside the horizon. We will now present the 
results for the standard and modified scenarios. All the relevant Whightman functions are listed in 
appendix IA.2I 



4.1.1 Standard scenario (review) 

In the standard scenario it is r] in = — oo, and the Whightman functions are obtained from the definition 
dU) using the mode functions in equation (|7|), see appendix IA.2I Inserting them in equation (j23j) . one 
finds @] 

(24) 

Adding the disconnected contribution in equation ()18p to obtain the bispectrum for £, and taking the 
squeezed limit fci <C %, &2,3 ~ finally 

(C^C^L ^ (2vr) 3 ,5(^ - n^fajP^s) . (25) 

i 

Equation (|25|) is the so-called local form of the bispectrum [3HB] . 



4.1.2 Modified initial vacuum 

As the computations are simpler in this case, we present it first. The time integration in equations 
(1191) . (|23p has now a lower limit at the finite time t]. = r] c . The relevant Whightman functions are 
listed in appendix IA.21 The result for the bispectrum for generic ki,k2,k$ has been first presented 
in [15] • Writing the result in the form (j2T)]) . (j2Tj) using ([231), yields 

aw, ^ = f: ^ = - e ^ + fc , 2 + t Re ( i - kh ~ kj)vc ^ w 

In the squeezing limit ki <C ks, &2,3 ~ the contribution proportional to f3 kl is very subdominant, 
as it is given by 

f*f * -Re[C*(l-^ sr ' c )]«l- (27) 

Larger contributions are found when the perturbations depending on the large momenta in the 
squeezing limit are in opposition of phase, that is, kj is either or k^. Indeed, in that case the 
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denominator and the argument of the exponential in 5F^ " are very small: using k^i.j = {k\ 
kj + 2k\kj cos 6 j) 2, from momentum conservation, and expanding in ^ = one obtains 



for j = 2 or 3 kh — kj ~ ki(l + cos#j) = k\ v e . (28) 

f— ^ fcl<fc2.3 1 

so that 

s 0=2,3) ^ _ k 1+ k 2 +h R r mv* , _ ^ ^ , 

1 fci«fc s fcif fl L J V /J 

It can be easily checked (for example by plotting) that the one in (|28p is an almost perfect approxi- 
mation for all values of 9j already for j~ = ^ ~ 0.1. 

The form of the correction <5.F 1 ' 7 ' 2 3 , then, depends on the interplay between the scales k\ and rj c . 
Recalling that we consider the realistic limit, where k\ is small but not zero, there are two asymptotic 
possibilities where the result can be evaluated most explicitly: 

1) \k lVc v gj | » 1 =* 5Ft 2 ' Z \ * -2^v 1 Re[)C] (30) 

3 fci<fcs ki 3 

2) |fcufc«,.|«l => 5F 1 (J=2 ' 3) ~ -2fe S 77 c Im[^r] ■ ( 31 ) 

3 fel<fcs 

We see that these contributions dominate over (|27p because of the large factors ^ an d ^s^c- In the 
case 1) we have eventually taken into account that in standard observables such as the halo bias, see 
section [6j the final result involves an integration/average over quantities such as 0j,ks and so the 
terms with large oscillations averages to zero. 

Thus, the leading correction to the squeezed limit of the three-point function for non-zero (3k reads 

tpiCcUCctfcijZ cub = A (h, ^ ^ S Fl ( kl ,k 2 , k 3 ) (2nf5(J2 h) £T cuh Pjh)P st (k s ) , 

(32) 

with 

3 

_mv, min cub „,mv, min cub 



^mv _ mv, mm cub mv, mm cub I ^ c fej "q. iLC L^fefi J \ n i.'IC"e j \ ^ * /ggN 

For observational purposes, for example concerning the halo bias, one is interested in the magnitude 
of these corrections and in their dependence on the probed large scale k^ 1 . We will discuss these 
aspects in details in section [5j We can already notice, however, that the /ci-dependence for the 
leading contributions (|32p . f|33[) is fully determined, since /3^J depends on the momentum ks and not 
on k\. The magnitude of the contribution will have instead to be estimated and bound using the 
phenomenological constraints reviewed in section [373] 



4.1.3 Modified dispersion relations 

The bispectrum in this case was first computed in [16], to which we address the reader for more 
details. The Whightman functions necessary to compute the three-point function are obtained from 
their definition in equation @ using the field solution (|10|) , Hence, the time integral in the correlator 
(|23p breaks up into the different intervals of validity of the piecewise solutions ()10p . In [16] , it was shown 
that the leading contribution occurs when the Whightman functions have support in the intervals IV 
and III, see figure [TJ Other contributions are indeed suppressed by A, defined in equation (jlip . or by 
a rapidly decaying integrancH, see [16]. 

5 The path of integration in time must be chosen such that the oscillating piece of the integrand becomes exponentially 
decreasing for ?y — > oo. This corresponds to taking the vacuum of the interacting theory [4]. 
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We concentrate, therefore, only on the leading contributions from intervals IV and III. It is con- 
venient to write the result in terms of the variables 

P^(V) _H X . = A_ (34) 

V- A - A k^V, x.-^. (34) 

In particular, when taking the squeezing limit, it is fc max ~ fc 2 ,3 ~ ks and xi <C 1. Since y is a rescaled 
time coordinate, sometimes we will call it simply "time" in this section. 

By inserting the relevant Whightman functions, listed in appendix \A.2\ in equation (|23p and 
writing the result in the form (|20|) . (|2ip . with (|24p . the corrections &F, 1 , read |16j 



i-^S'o({a: h ^ J },2: j ,j/ / ,m) 



(35) 



5i? m ={o,i}( x i' x 2,a;3,y)=X) le (^T ^17 / d y' x t9{{ x h^j},Xj,y',m)e l H 

j=i L 7 % 

where x t = X^=i x «> an d 9{{ x h^j},Xj,y' ,m) is reported in equation ()106p of appendix I A. 31 So is 
defined as^| 

S ({x h ^},x j ,y',m) = f V dy"(Y,u(xh,y") + (-l) m u(x 3 ,y")Y h,j G {1, 2, 3} , (36) 

The limits of integration for the variable y' in (|35p have been discussed in [16], and we refer the 
reader to that reference for the details, reporting here the results. The upper limit is y ~ because 
the bispectrum is evaluated at late time r] ~ 0, while the lower limit y n is determined by the smallest 
among ??i fcl=1 ' 2 ' 3 \ which bound region III for the three momenta k i=1 23 , see equation ()10p . As shown 
in |16j . it is y n = ■^■^max r ?i fcmax ^ ~ — 1- This is indeed straightforward: the corrections to the linear 
dispersion relation at the time rj^ that separates region III and II must be quite important in order 
to drive the frequency close to the turning point (see figure [1]), and so, since these corrections grow 
as powers of see ([9]), it must be |j?(7y^)A -1 | ~ 1. Furthermore, since pfaW) = -krj^H, the 

smallest boundary time rj^' is obtained for the largest k. Thus, we have to consider r/( fcmax ) and, from 
the definition of y, it must be |yj = |p m ax(??j; fcmax ' 1 ) A" 1 ) ~ 1. Observe then that \/\y\ < \y tt \ it is 
< I x i2/J 1 in the squeezing limit, so effectively F(^kir)) = F(x\y) ~ 1 and uj(x\y) ~ x±. 

The integral is then divided in two contributions, according to whether the Whightman functions 
depending on /c2,3 have support in their region IV or III. In the squeezing limit, they will have support 
mostly in region III, because the modes will become superhorizon (region IV) only at rjs — —k^\ — 
~ tl> that is for y$ = — ^ ~ 0: this time is also pretty close to the late time y ~ that the 
bispectrum is evaluated at, thus the contribution from when all Whightman functions have support 
in region IV is subdominant. We concentrate now on the contribution from region III (of course ^ 
will have supports both in its III and IV regions). 

The one in equation (|35|) is a typical Fourier integral, dominated by the contributions of points 
where the strong suppression (-^ 3> 1) from the oscillating phase of the integrand is reduced. Con- 
cerning SFq, the only such points are the boundary ones: in fact, there are no stationary points since 
the first derivative of the phase So (£1,2,3 j jA m — 0) is the sum of the positive-defined frequencies and 
thus is always positive and never zero, see equation (|36p for m = 0. Fourier integrals like this can 
be approximated for jj 3> 1 by integrating by parts (see [27] and section 4.1.2 in [IS]). Thus, in the 
squeezing limit x\ <C 1, x 2 3 ~ x s — 1 we find 

SF ,* 1- J Im[g(x 1 ,x 2 ,x 3 ,-l,0)e^ So ^^'- 1 '°)]| Xl ^ ~1, (37) 



(a;i<l) 



u 2,3~ 



J Observe that uj(x,y) = xF(xy) is dimensionless. 
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where in the last passage we have taken into account, as before, the averaging to zero of the rapid 
oscillations in observables such as the halo bias. As we see, this contribution does not lead to very 
significant corrections to the standard local form of the bispectrum. 

We expect a more interesting behaviour for what concerns 8F\ |16j . In the squeezing limit x\ <S 1, 
we find that the dominant contributions are those where the mode function in the negative-energy 
branch is the one depending on ki or ^3, because in those cases the suppression from the oscillations 
is very reduced, since the perturbations depending on the large momenta are in opposition of phase: 



for j = 2 or 3: 



So({x h ^j},Xj,y',l) = I dy"\y2,u(x h ,y")-u(xj ,y") J ~ xi j dy"(l + d x u(x, y")\ cos Oj) +0{x\) 



Xl y'(l + F(y')) cos(%) + 0{x\) = xi v e {y') 



(38) 



where we have used equation ^ and x s ~ 1. 

The solution for the integral has two asymptotic regimes, depending on the magnitude of the phase 
factor A x\ and the behaviour of the integrand at boundary and stationary points, see [27] , We have 
listed in appendix IA.31 table [H the relevant information, which leads to the solution 



6F, 



3 Re 



■-E 



J=2 



l+cos(6».,) 



A^— 4r(^) 



1 \ Im 



1 



mdr* i 



* . n signal 



1 K+l 



(k+1)|F< k )|^FT 
3 



>x- 



11 u B j ^ A li 



3 Im 



--£ ■ 

3=2 



C 9{-De 



► if^x^l 



+ E E(£) "* 

j=2 {y,} x ? 



2Al mE [^(0(l) +0 (x 1 ))] 



J'=2 



if ^S!<1 



This complete leading- asymptotic formula looks complicated, and we write it only this time for illus- 
trative purposes, but the final result will be simpler and clearer. The notation is explained in appendix 
IA.31 Let us recall here the main points. We have defined 



3 TVM±. 

z£*' v *Xv*) = E 



(39) 



where, v[ u *\y*), y **" 1 '^*) are the coefficient of the leading behaviour of the functions v g , (y), g({x},y) 
in the squeezing limit in proximity of a stationary point y*, see table [TJ The sign ± in the phases 
depends on v e {y) — C„.(y*) being even or odd under (y — y*) — >• — (y — y*), for y ~ y*- If ^* is an 

integer, = n! -1 <9™,?) e . (similarly for //*). Recall also that u e . = l+cos(0j). 

Finally, the contribution of the (nearly) folded configurations (singled out by the condition v g . = 
1 + cos(#j) < is determined by the lowest order correction to the dispersion relation, which, 

expanding ©, would go (with model-dependent ~ 0(1) and re) as 



A' 



(40) 
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In equation (|39p . the contributions in the first line for jjX\ 3> 1 arise from the boundary points, 
the second line from the possible presence of a number of stationary points {y*}, which for generality 
we have taken as stationary of order z/* for v e .(y), and as zeroes of the function g of order see [27] 
(/i* and Vif are generally different for the various points). The presence of stationary points is strictly 
dependent on the specific dispersion relations, that is the specific model. However, their contributions 
will be suppressed in the final observables such as the halo bias, because of the averaging to zero of 
their large oscillations, visible in (|39p . 



Inserting these results in equation (|20p , and using (|2ip , (|24p , the leading correction to the standard 
consistency relation due to non-zero fik reads 



(41) 



r b = e 

mm cub * — 

i=2 



Im 



Q mdr* ^BignCFW) 



(k+1)|FC=)|^TT 



ifl + cos(^)>f i 



ifl+cos(^)<f^ 



if 4sci » 1 



if Axi < 1 

(42) 

where we have accounted for the averaging to zero of the rapid oscillations in (|39p in observables such 
as the halo bias. Again, at leading order the dependence on k\ of the result is fully specified, and 
there appear enhancement factors such as -jj and ^ • We will comment on the result in section [SJ 

The final result is very similar to the one obtained in the modified vacuum case, which points 
at the fact that indeed the squeezing limit is dominated, at leading order, by the generic features 
common to both scenarios (such as particle creation, interference effects). The only difference is that 
in the modified vacuum case, the folded configurations yield a contribution corresponding to k — > oo 
in ()42p . since the phase of the integrand and all its derivatives are identically zero, see (|33p . In the 
case of modified dispersion relations, instead, the true analogous of the folded configuration (which 
would be ^2h^j u ( x h , y') — uj(xj , y') = for some x's but for all y' 's) does not occm|j and k is finite. 



4.2 Example 2: Quartic derivative interaction 

In the standard scenario, the contributions to the bispectrum from higher derivative interactions are 
suppressed in the squeezing limit by powers of [B]. It is interesting, then, to study if the additional 
contributions to non-Gaussianities due to particle creation in the scenarios of modified vacuum and 
modified dispersion relations would lead to the same suppressed contributions to the bispectrum and 
the same overall scale dependence. 

In this section, we present a detailed example of the kind of calculations involved, before passing to 
the more general discussion in section I3~3l We take as our representative example the cubic interaction 
coming from the correction £ HDI = V - detGg^((V<!>) 2 ) 2 to the effective action of the inflaton $, 
where G is the metric, g the coupling. In the standard single-field slow-roll scenario, [21] has discussed 
that this is the lowest-dimensional operator in the class of operators which is the most important for 
non-Gaussianities. Nonetheless, for the standard scenario its contribution is very suppressed in the 
squeezing limit, as we have said. 

7 It can be easily verified at least if ui is given by elementary functions and/or power series. 
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Expanding in perturbations and converting in terms of £ as explained in details in |15|.ll6j . from 



£ HDI one obtains the cubic term 

%) = " / d 3 xa^^C(C 2 - (9 iCc ) 2 ) • (43) 

Because of its origins from the term £ HDI , which contains quartic power of derivatives, we call this cubic 
interaction "quartic derivative interaction" . It is allowed to treat these interactions perturbatively also 
when the physical momenta p = a~ 1 k are such that ^ is not small any more, because, compared to 
the operators changing the dispersion relation, they are further suppressed by additional powers of 
the perturbation fields. 

The bispectrum, from equation (fT9j) . reads in this case 

(C,fc 1 (^C CiA : 2 (r/)C Ci fc3^))= 2Re[-i(2vr) 3 ^ 3 )(^^)^ I /' r ' d V ' a{d v G kl (r,, V ')d v G^( V , V ')d v G k3 (,r,,r,')+ 

+ h ■ k 2 G kl (77, rj')Gk 2 (??, r/)d n G k3 (77, t/)) + permutations , (44) 
where, again, 77 ~ is a very late time when all modes ki are outside the horizon. 
4.2.1 Modified initial state 

In this scenario the leading contribution to the bispectrum for generic external momenta has been 
computed in [9j[15]- Writing once more the result in the form (j20|) . (f2Tj) . with (f24"|) . it is found 

, (45) 



5F 1 (k 1 ,k 2 ,k 3 ,i] ~ 0) = Re 
with 



ViMa E C f dv'e^ Pjfa k 2 , k 3 , V ') + c.c. 

ft l fc 2 fc 3 j—i Jr)c 



Sr>i> ^r 2 "? 4A 



3 ( 

Pj(ki,k 2 , k 3 , rf) = -kt Y[(^2 kh ~ k ^ ~ i7] ' ( ~ fc i ~ ^ 2 +2) fc j+i( fc i+2 - + (^+2 - ~ tf+i) 

i=l /i^i ^ 

\ 3 
k] +2 (k j+1 -k 3 )+(k]-k] +1 -k] +2 )k](k 3+1 +k J+2 ) J +(-r/) 2 (^ fcfc-fcjOCjI h){k1-Ak j+1 k j+2 ) , 

(47) 

where j is defined modulo 3. The factor C has been written in terms of the variables Xi defined in 
equation (f34"j) . 

As in section 14.1.21 we find that in the squeezing limit the dominating contributions to the integral 
occur for j = 2, 3 in the sum in equation (|45p . in particular, the largest ones come from the terms of 



order rj and rj 2 in Pj. Their integration yields 



, yv t lw, l + e^'V-l + iJfeif/cW- ) j« 
d,/ e *(£h* fc fc -* 3 -)»J ,/ = — i- = A- , (48) 



k\v 2 g k\v 2 ' 



/•o 

/ dtfe % 

Jr\ c 



(2w*j k h -k 3 )r, 12 = . , j_ _ 2. 4g 
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so that 



§F ~ > — Planck^. -2 

J=2 



M/Ci? ]-co S (^-)Re[C7 : 



* -0) 
1 



(50) 



Thus, looking at (|48|) . (|49|) . we find again two asymptotic behaviours for the correction to the 
squeezed limit of the three-point function depending on whether \kii] c v g . | 3> 1 or \k\r] c v e | <C 1: 



quart dor 



(2vf6(J2ki) B™ ItdcI P st (h)P st (ks), 



(51) 



with 



Z3 



quart der 



gH 2 M 2 

g Z Planck 



2c 



gg2 M Planck fcl rt-^-n ^2 



if |fci77 c u e . | > 1 



F-^fe) 2 E I fe(*fiifcK, Im tC ] +«»&) ] if \k lVc v gj | « 1 . 

j=2 v 7 

(52) 



In the first of (|52|) . we have again taken into account the averaging of large oscillations to zero in 
observables such as the halo bias. Recall that |/c5ry c | 3> 1 (in the NPHS scenario, in particular, 
\ksVc\ = 37 )■ We comment on the actual magnitude and scale dependence of these corrections to the 
bispectrum in section [5) 

4.2.2 Modified dispersion relations 

The three-point scalar correlator for generic k\ , /C2 , &3 in this scenario was computed in [16] . We focus 
on the most dominant contribution, which occurs when the Whightman functions have support in 
regions IV/III, see (|10p . as discussed in [16] and reviewed in section 14.1.31 Using in equation (|44p the 



Whightman functions listed in appendix I A. 21 and writing the result in the form (|20|) . (|2ip with the 
variables in (|34|) . the leading correction linear in |/3^ dr | is again given by a Fourier integral, see |16| : 

/A x3 



6Ft 



dy'y' 2 e i £ s o({nh}r>:j,v'A) 



<({x h },Xj,y'] 



(53) 



where the limits of integration are those presented in section 14.1.31 The function q is reported in 
equation (jll)9|) of appendix IA. 31 while So and C are defined respectively in equation ([55]) and ([4*6"|) . 

When taking the squeezing limit, the dominant contributions are again those where the pertur- 
bations depending on A;2,A:3 are in opposition of phase. Then, $o ({^1,2(3) }> ^3(2) jZ/i 1) behaves as in 
equation ([38]) . The integral can be solved by asymptotic techniques, and we obtain at leading order 

6 ^C,kSc,kAkXuL der k ^ ks ( 2 ^) 35 E^Cl dor P nt ( k l) P st (k S ) 



with (recalling appendix IA. 31 and table [Q) 



B 



mdr 



gH 2 M 2 

_A(:Z- Planck 

A 4 



quart der 



3 

xE < 

J=2 



fc^ r(3)+cos(6> j ) 
fcl (l+cos(6>j)) 2 



Re[/3 



mdr* 



H\Hk s ) (k+1) 



(#) 2 ("fe 0(1)] + cos(0,)Re[C* 



(54) 



if ^xi > 1 



if ^xi < 1 
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where, again, we have taken into account the averaging to zero of the strongly oscillating terms in the 
final observables as the halo bias. Equation (|54p matches ([52]); in particular, note that in the modified 
vacuum scenario, as we said, the folded configuration has k — > oo, and the dominating term in the 
second line of (|52p is the one proportional to Re[/3j^ v ], since \k\r] c v e | <C 1. 

4.3 General analysis 

Let us first summarize the key-points concerning the violation of the standard result for non-Gaussianities 
in the squeezing limit, which have emerged from the detailed examples discussed above: 

• non-Gaussianities can be enhanced^] in the limit by two effects: 

— interference, which reduces the suppression due to the oscillating phase of the integrand 
in the time integral in the formula of bispectrum, and is strongest when the perturbations 
depending on the large modes ^2,3 ~ ks are in opposition of phase, 

— accumulation in time, which leads to larger enhancement from the time integration when 
interactions scale with higher powers of ~ 7/ 

• if the probed squeezed mode and the time/scale of new physics are such that \kir] c \ <C 1, 
or ^77 <C 1, then the perturbation depending on k\ is outside the horizon already at the 
time/scale of the new physics, and thus the result is still of the local form but with enhanced 
non-Gaussianities (because of the interference and particle production for the other two modes 

*2.3) 

• if instead \k\r} c \ > 1 or 7^:77 > 1 (in a realistic squeezing limit k\ 7^ 0), the perturbation 
depending on k\ is inside the horizon at the time/scale of the new physics and the bispectrum 
is not of the local form 

We want now to generalize our analysis beyond the two examples of interactions discussed above. 
It is convenient to adopt a different gauge than the one we have used so far: in the new gauge it is 
C = 0, and the perturbations are accounted for by a "matter field" perturbation ip. This gauge is fully 
equivalent, of course, to the one used before, and it is useful because we will have to consider only 
the matter action when expanding in perturbations. In particular, we can choose ip to be the inflaton 
perturbation, expanding the inflaton in background plus perturbations^! 

$(x,t) = (l>(t) + <p(x,t). (55) 

The results in terms of the variable ( are then obtained at the end by performing a gauge transfor- 
mation as described in [3]: 

C = -H% + 0{e 2 ,cp 2 ). (56) 

9 

so that 

H 3 r 1 

(CCC) = — —( t P l P l p)+ (C(e 2 , </9 2 ) (^99) + permutations . (57) 



The actual magnitude of the enhancements depends on the precise value of \j3k s | and the scales of "new" physics, 
and thus can only be estimated, as we will show, in our phenomenological approach. 

9 This is quite intuitive, since the contribution we are interested in will arise at early times (large \r)'\) around r\ c for 
the modified vacuum scenario, or around the time of WKB violation/particle creation 7711 ~ Ak^ 1 !!^ 1 for the modified 
dispersion scenario. More specifically, in the cases where |fei?7 c | -C 1 or fj-jj -C 1, since the integral for the bispectrum 

becomes approximately of the form drj' rj' n ~ r]"^, one can easily see the growth with n. We will soon show a similar 

growth of the enhancement also in the cases where \kirj c \ S> 1, jr"]j 3> 1- 

10 This gauge also makes clear from the outset the actual order in slow-roll parameters of the cubic couplings, see [4]. 
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In particular, it is enough to consider the gauge-transformation at leading linear order, because the 
higher-order corrections will not yield time-integrated contributions to the bispectrum at tree-level, 
and, as it appears from the previous sections, it is the time-integrated contribution that leads to the 
leading modified results (see also [9][TI^[T6] ) . 

The most generic effective action for the inflaton <I> in single-field slow-roll inflation has been 
presented in [28] 0. In a Lorentz invariant theory, the terms in the action will be function of the 
Lorentz invariant objects <£, g^ v d^d u ^, Wd^d^, R. In particular, those that are functions 
of $ only, and not of the objects involving derivative, constitute the scalar potential V(3>). Recall 
also that all second order time derivatives (and first order derivatives of auxiliary fields) in higher 
order terms in the action (more than quadratic in the field) should be eliminated using the equation 
of motion, see [2H]. This applies for example to □<£ and W" . In a Lorentz-broken theory (such as it 
is the case with modified dispersion relations), there are also additional terms that are not built up 
only from the Lorentz-invariant operators listed above. These terms involve higher derivatives of the 
fields that can be written as purely spatial derivatives in a convenient system of coordinates o 

By expanding the effective action in the inflaton perturbation ([55]), given the fundamental building 
blocks listed above (including the higher-derivative Lorentz-breaking ones), and using conformal time, 
the cubic couplings for the perturbation ip have the most general and schematic form 



/ 



drjSx a 4 - = — T- 58 



where d can be either a spatial or time derivative (comoving coordinates) recalling however that in 
the case of time derivatives, their order is at most 1 (as we said higher orders must be eliminated by 
using the equations of motions, see [28]). We have sloppily indicated both time and space derivatives 
with the same symbol here for notational simplicity and because it will simplify the counting of powers 
of momenta in the correlators, but we will distinguish more carefully the two kinds of derivatives in 
the calculations. Finally, the dimensionless coefficients A nms can depend on <j) and in general on the 
slow-roll parameters, and are constrained by backreaction. 

Let us comment on the possible values of n,m,s in the cubic couplings. In a Lorentz-invariant 
theory, given the Lorentz-invariant objects listed above, n, m, s can be at most equal to 1 0, while 
in a Lorentz-broken theory their values can be higher (for the higher spatial derivatives couplings, 
in the convenient frame). Observe that in a gauge-fixed theory, when solving for the gravitational 
constraints (for example those related to the lapse N and shift-vector N l in the ADM formalism) there 
will also appear non-local operators. In the gauge we are using, they arise because of the solution^! 
N i ~ ^(d)" 2 ^ + 0(e 2 ), see @|. This means that, when expanding in perturbations, n,m,s can 
also assume value -1, because of the combination (d)~ 2 <p, in the leading order contribution from N l , 
relevant for the bispectrum. However this combination is acted upon by d z , which compensates the -1 
with the additional derivative. 

Observe also that N % enters the action only via the metric g^ u or the extrinsic curvature. In the 
first case, the spatial index of N l will result always contracted with a spatial derivative, as it can 



11 Another way of writing the most generic effective action for a "matter scalar" takes advantage of the so-called 
Stueckelberg trick using the Goldstone mode for the broken time diffeomorphisms in the unitary gauge, see [29] . 

12 Lorentz-breaking always implies a privileged frame, the convenient system of coordinates is precisely the one adapted 
to this frame. 

13 Recall that, since we are considering only scalar perturbations, after the expansion spatial derivatives indexes must 
be contracted by the metric, leading thus to contributions proportional, by permutation, to ki ■ kz, fei • k3, hi • fe. 
In particular, in the squeezing limit the contributions proportional to fci ■ ki and k\ ■ k$ give an overall contribution 
k\ ■ (&2 + £3) ~ k\ [BJ, which is subdominant with respect to the contribution ki ■ hi that goes as ki ■ k-$ — > —k2ks = —k% 
in the limit. The latter is thus the one we consider as leading in the following. 

14 Obviously, integrating by parts, derivatives can be "moved around", in which case the more appropriate condition 
is that max(n + m + s) — 3. What we will say is valid, of course, also for the terms obtained by integration by parts. 

15 The standard notation (d)~ 2 for the inverse of y^-d^dj is formal. It becomes clear after a Fourier transform. 
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be seen from g^ u d^dy^ considering that the shift vector enters the components g° l = N~ 2 N % and 
gV = N^WNi . Also when coming from the extrinsic curvature, a component N l will always be 
accompanied by a spatial derivative, because the extrinsic curvature is defined as = (2N)~ X (juj — 
2 (3) V {i N j} ) = (2N)~ l {hij - h ik djN k - h jk diN k - N k d k hij), where h ih (3) V are the ADM three-metric 
and relative covariant derivative. These properties imply that in a coupling where originally N l was 
present, the combination d l (d)~ 2 (p will multiply at least one other factor with a spatial derivative. 
Furthermore, couplings with n = m = s = arise from the expansion of terms such as 

V"'(6) o /<9 r $\" , o 

V($) ^ —^-tp 3 , or (— ) d> m -> (dl<t>) u r- V , (59) 

and thus are higher-order in slow-roll and would give a negligible contribution to the bispectrum (see 
also [4"ll28tl2"9"]). We will neglect them in the following. Moreover, couplings of the form (p(d(p) 2 , 

where d denotes spatial derivatives, are reabsorbed (eliminated) by field redefinitions such as <j) — > 
4> + ^j^c/) 2 . Finally, couplings with a single time or spatial derivative are absent because of isotropy. 

The two detailed examples that we have analyzed previously fall into the scheme we are describing. 
In particular, the "quartic derivative" coupling (|43p derives precisely from the inflaton action, as 
discussed in section 14.21 and is actually the sum of two elementary cubic couplings (one with only time 
derivatives, the other with one time and two space derivatives). In the "minimal cubic" coupling (|22p . 
instead, we find the combination {d)~ 2 d r} that we have been discussing. Summarizing, once written 
in terms of the inflaton perturbation, those couplings are of the type 

. r -ir- i min cub 

minimal coupling : {n,m, s\ = {1,1, A n 1 = v2e— (60) 

Planck 

, , quart dcr fe HM m , , , 

quartic derivative : n = m = s = 1, A ni = W - g ^ 2 anc . (61) 

Finally, not writing the momenta conserving delta function, 2tt factors and overall minus signs, 
the contribution to the bispectrum from equation ([58]) for a generic cubic coupling reads 

(Cfc^C^h K ^ A n+mT B -i f dr l' <Vr~ n ~ m ~ s (d n G) kl (d m G) k2 (d s G) k3 + permutations 

+ complex conjugate (62) 

where we have used equation (|56p to cast the result in terms of the perturbation £ (the Whightman 
function in this formula are those of the (^-variable). (d r G) k indicates the Fourier transform of d r G. 

We now study the dominant contributions in the different modified scenarios, recalling that, as it 
has emerged from the detailed examples, it occurs when the interference effects are stronger -that is, 
when the perturbations depending on the large modes ^2,3 ~ k$ are in opposition of phase-, and for 
the terms in the integrand with higher powers of r( (higher powers of - ~ r/) and, obviously, lower 
powers of k\. We will not pay attention to numerical factors of order 0(1) in the formulas as they are 
irrelevant for our considerations. 

4.3.1 Modified vacuum, Lorentz unbroken 

The relevant Whightman functions go as, see appendix IA. 21 

tr4 ±ikT}' ttA p ±ikr)' 

(fyG^MO,*/) = ±-p-fcV^T> (^°G ± ) fe (0,r ? = ± (± iki Y(l T ik V ')^-, (63) 
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so that the leading contribution to (|62p . according to our analysis, is 



P \^k\^k c 2 ks ' A A n+m + s_1 ^ 2 J 0323^3^3^,3 a ° J v ^ 



gifci (1+cos 0j )r;' 

~ h hermitian, (64) 

2 n t 



where we have used equation (|28p and a,b,c run over 1,2,3 accounting for the permutations, from 
which we take the dominant ones according to footnote [TBI Finally, nf = {0, 1} accounts for the 
presence or not of time derivatives acting on G(ki,rj'). 

We concentrate now on the time-integral, which can be computed in closed form: by writing 
v = n+m+s(+nf — 1) and recalling our discussion on the values attainable by n,m,s in section 
we find 

r/r/ r/- x +° pifciCl+oMfliV = (ST ( V ~ l \ r ^- ik ^) V ~ l ~ r Jk^Al+cose,) _ 

„ c V V ~ {ikl) v{2^{ r ) (i + cos e .y+i (1 + cos fly)" 

f-f if ^(l+cos^Kl 

~ I W if|Mc(l+cos^|»l. 

where in the last passage we have written the leading contribution, also considering the averaging to 
zero of large oscillations in the final observables. 

Inserting this result in equation ([64"]) . we obtain the leading contribution 



hikkCtfT™ *T pmpm (66) 

with, having substituted back n + m + s in place of v, 



B 



Pks V2~e H 2 y k 



h\ ,l+min(n,ra,s) 

5 



H%\k S ric\y 



(C - (1) + rf^ ) if| ^'« 1 



3 

+ complex conjugate (67) 



S(-»fxfP" +s (;M^+ :?( J if|Mc^ l»i 

r — ' x iV fi/L / x U/i n + m + s + n: -1 ' 3 

J=2 »j „ * 



where 



u 9j = (1 + cos c nma (x) = 2x[n+ ^l s _ 1+x) , d nms {x) = -{-\) x {n+m+s-2 + x)\d x , (68) 

and c\,d\ are order one numerical factors, possibly depending on 8j. 

We observe that the generic contributions present both enhancement factors (such as |fcs?? c | S> 1, 
or I s - S> 1) as well as suppressing ones (such as (#) n+m+s anc j A nms ). We also observe that for 
\k\rj c \ S> 1 the contribution will not be in general of the local form. However, such contributions could 
be suppressed and thus subleading. We will discuss enhancements and scale dependence of the results 
in details in section [5] Observe that for the couplings ([60]) . ([6"T]) . equation (f6T|) leads to the results we 
have found before, with the proper c\,d\, see equations (f3"3|) . ([52]) . 

16 Recall that we are not considering the terms with n = m — s = because they are subdominant in slow-roll, that 
couplings with only two spatial derivatives are reabsorbed by field redefinitions, while those with only one derivative of 
any type are not present because of isotropy, and also that in the case of time derivatives their order is at most 1 (higher 
orders must be eliminated using the equations of motions, see [28]). 
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4.3.2 Modified dispersion relations, Lorentz broken 

As discussed in section 14.1,31 the leading contribution to the bispectrum arises when the Whightman 
functions have support mostly in the interval of times in regions IV/III, given the solution (jlOp . The 



Whightman functions go as, see appendix I A. 2 1 



(^G ± ) fc (0,r / ')=-i T -fe 



H 3 ftM)^* ^) 



a(r]')2k 2 



H 3 yWfA; ri) e ± ^ n(fc ' r?,) 



a{rf) 



where X i^iV') is defined in equation (|103|) . 

Once again, it is convenient to write the contribution to the bispectrum in terms of the variables 
defined in equation (|34p . The leading contribution then reads 



A, 



J A3 



J=2 

dy' {-y') n+m+s - 1 (i7 ( *W 



2 3 k 3 k 3 k 3 \H 



in+m+s 



« e 



+ hermitian, (70) 



where we have used equation (|38p . n[ a 6 c) = {0, 1} accounts for the presence or not of time derivatives 
in the coupling acting on Ck abc i an d a,b,c run over 1,2,3 taking into account the permutations, 
from which we take the dominant ones according to footnote (|13p . Now, depending on the specific 
Lorentz-broken model, n, m, s can assume values larger than 1 (n + m + s larger than 3 1^1). 

The integral can be again obtained by asymptotic techniques so that the leading contribution reads 



where, having written n + m + s = v and used table [1] in appendix IA.31 



A„ m „ AM D ,„„„,, , k\ v l+min(n, m,s) 



3 

i=2 



Planck / 1 \ J 

.27 # 2 { k s ) 

(§r"(o(i) + (#fe) nf " 1 o(i)) 



i« sign(F( K) ) 

/flt. s ^ r ^) eT (K+1) ri , r M ( Hks< 



if AAl i 



if AiiL -*> 1 
11 H k s 



> + 



(71) 
(72) 

complex 
conjugate 



where 



■ 7r(K+r+n) signffjj) 

tt i 00 ui n/ ii+K+f+n u' 2 

1 r,n=0 1 



r(_^)|F(^|^- 



(73) 



contains all the subleading asymptotic contributions (the c r>n are calculable order one coefficients, 
whose precise form is irrelevant for us, obtained from the expansion of the integrand in (|70p ). 

This result is very similar to that obtained in the modified vacuum scenario, except that now n, m, s 
can also assume values larger than 1 in the case of couplings deriving from the Lorentz-breaking terms 
in the action, and that k — > oo for the folded configuration in the modified vacuum scenario, whereas 
here, for the case of modified dispersion relations, the contribution is again determined at leading 
order by the first correction to the standard dispersion relations (see (|40p ). A part from this aspect, 



However, time derivatives are at most of order 1, see footnote 1161 
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the same comment as the end of the previous section applies. The similar structure of the two results 
can be seen as showing that indeed the squeezing limit is dominated, at leading order, by the generic 
features common to both scenarios (such as particle creation, interference effects). 

For the couplings ([60]) . ([6T]) . equation ([72"]) leads to the results we have found before, see equations 
(|42p. (|54p . Note, in particular, that since the one in equation (|43f) is the sum of two elementary cubic 
couplings (one with only time derivatives, the other with one time and two space derivatives) such 
that for the folded configurations the leading contribution is cancelled out, in those cases the result 
([31]) is given by the first subleading correction in equation ([72]) . 



5 Signatures of very high energy physics in the squeezed limit 

We will now comment and discuss the results found so far. For observational purposes (regarding the 
halo bias, for example), the most interesting features of the squeezing limit of the bispectrum are its 
dependence on the large probed scale k^ 1 and the magnitude of the non-Gaussianities [7J[8]: therefore, 
we will focus on those aspects. 

From our analysis, it appears that the squeezed limit of the bispectrum shows different kinds 
of behaviour regarding /ci-dependence and amplitude, subject to the interplay between the squeezed 
momentum scale k\ and the scale/time A/rj c of new physics. In particular, the different behaviours 
occur depending on whether \kir] c \ or f|-7j- are larger or smaller that 1. 

We will momentarily investigate what the different modified scenarios predict concerning the value 
of | A;i?7 c | or t^jj, but let us postpone this question for a moment, and discuss first in more details what 
(and really if) significant enhancements could appear in any case, and what would be the dependence 
on k\ for the bispectrum in the squeezing limit. 



5.1 Enhancements 

We start by discussing the enhancements. Since we have adopted a phenomenological approach we do 
not have the full-detailed knowledge of the magnitude of j3k s , which the leading contributions to the 
bispectrum depend upon. However, we know the constraints that the Bogoljubov coefficient has to 
satisfy for phenomenological reasons and self-consistency of the theory (see section 13. 4p , and thus we 
can estimate the largest possible enhancements. Observe also that the leading contribution depends 
only on (3k s , which does depend on the scale kg, but not on k\, thus the fci-dependence of the leading 
contribution is not affected by our ignorance of 

Cases \k\T) c \ < 1 and |^ < 1 

In this case the correction to the bispectrum has a parameteJ^l Jml ~ £>™ v ' mdr , where 



|B m - md 1 x ~ £ |AJ^^| r ^(|l) I+m1 "'"''"-' (^)" +m+ '(l + (^)"™- 1 0(l)) (74) 

< .via ^ %(|)— » (1,(^)^(1)) . m 

In the inequality (j75j) we have used the constraint (]16p on |/3/% g |. Here, D = \ksr) r \ in the modified 
vacuum case, while D = jj in the scenario with modified dispersion relational. The similarity of 
the results in the two different scenarios is an indication that the squeezing limit is dominated by the 
generic features common to both of them (particle creation, interference effects) . 



18 As before, 'mv' indicates the modified vacuum scenario, while 'mdr' the one with modified dispersion relations. 
19 However, in the NPHS approach |fcs»7c| = tj- 
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If we consider first D = \k$r] c \, the inequality reads 

IB^^'L < ^>/i^T (^i-H-iK^™..) (JP S ^E^)"^^( 1+ (J^^)-«°- 1 C p (1 ^ , (76) 

where p sw (r] c ) = — k S{1) rj c H is the physical momentum associated with the wavenumber at the 

boundary time r\ c . 

When instead D = such as for the cases of NPHS vacuum or modified dispersion relations, the 
expression for |i3|^ PHS ' mdr simplifies and we obtain 

NPHS,mdr r— A^ , h , l +m in(n,m,a) / ( h A. \ «f "1 \ 

1 L " VW V2 HA [ k s ) \ \k s HJ U[L) J- {N) 

From ()76p . ([77]), we can understand various characteristics of the amplitude of the non-Gaussianities 
when |/cir/ c |) <C 1. First of all, we see from (|76p that the maximum amplitude for modified vac- 
uum scenarios occurs when the physical momentum p s (j] c ) associated with k s at the particle creation 
time r] c is really at the largest scale that we can trust our effective theory at: the scale of new physics 
A. But this is precisely what occurs in the NPHS scenario of vacuum modification, leading to (|77|) . 
In the case of modified dispersion relations, equation (|77p arises for the same reason: the physical 
momentum that would appear in place of p s (rj c ) in the analogous of (fTB"j) is p s (r]^ s) ), at the particle 
creation/WKB breaking time r]^ s , which has in fact a magnitude of the order of A. 

However, in order to see if there can be actual overall enhancements, potentially interesting for 
observations, one needs to take into account also the other factors in the equations for Z3™ v ' mdr : both 

the suppressing (A nms e-y/|/I[ 1+min ( n ' m ' s )^ anc j ^he enhancing ( — ^ lck ) . In this respect, it is imme- 
diately evident that terms with higher number of derivatives are more suppressed (since min(n, to, s) 
will be higher). This is particularly true for the Lorentz-breaking couplings in the scenarios with 
modified dispersion relations, where n, m, s can assume values even much higher that one. 

As it was in fact imaginable, the cubic couplings (|22p and (|43p that we have discussed as detailed 
examples, which have the coupling constants and order of derivatives as summarized in equations 
([60]) . (f6"T|) . provide (among) the least suppressed contributions, with (|^-) and (^-) suppression 

respectively (a part from the suppressing factors A-^ , A m , of course). The minimal cubic 

coupling is in fact also leading in the standard scenario, see [6]. 

If we concentrate on these dominant (less suppressed) couplings, which we have studied in details, 
we find that there can be actual overall enhancement of the non-Gaussianities, with respect to the 
standard scenario, for certain values of the scale of new physics. Let us discuss a concrete example: 
take e ~ ~ 1(T 2 , H ~ l(T 5 Af planck , compatible with the results of WMAP, and A ~ l(T 3 Af planck , 
which corresponds to the super symmetric GUT scale. Then, if we can probe down to < 10~ 2 , 
the largest possible enhanced amplitude with respect to the standard scenario that we could observe 
would be roughly of the order of 

i „NPHS, mdr , , . 

\B Ucub~10e ~10(l-n s ). (78) 

in the case of the minimal coupling (|60p . while 

■ „NPHS, mdr , □ 

\B Ur t der ~ ^9 (1 " ~ 10 (l - U s ) (79) 

ion ()6ip with g ~ 10~ 2 . It is notable how, for the 

1 ' \ : i ] ] ! ' k 



in the case of the quartic derivative interaction (|6ip with g ~ 10 . It is notable how, for the higher 
(quartic)-derivative interaction, the additional suppressing factor (^") 2 ; and enhancing factor — -jp 
(essentially due to the power of a -1 scaling) compensate each other to give a result of the same 
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magnitude as that of the minimal cubic coupling: a result which is quite different from the one in the 
standard scenario (where the a -1 scaling of the interactions does not have any important effect, since 
the dominant contributions to non-Gaussianities occur at late times). 

Observe that, according to the analysis of [6], values such as ^ ~ 10~ 2 are indeed within the 
reach of LSS investigations, such as for example EUCLID, which in the most optimistic estimates can 
probe values down to ^ < 10~ 3 [6]. 



Cases |ki77 c | > 1 and » 1 



In this case the correction to the bispectrum has a parameter /nl ~ B™' mdr with 



\B 



~ e 



I a I ^nms ^jjglanck / ^1 \ l+min(n,m,s) „( n m s ) ,H ks, 
\ lma ^Planck / ^1 -\l+min(n,m,s)p( nimi g) % 

V2 #A l fe s ; l AV 



(80) 
(81) 



where in the last inequality we have used the constraint (|16|) on |/3fc s |. Here, 



E (n,m.,s) r^_^S_\ 

A fci 



A fci 



> 



(^) n+m+S (l+(^) n *'-b(l)) 

(f M ) ^ (l+E oo =o0((f Mp )) mdr _ 



_H_ fcs 

A fci 



if 1 +008(0,) < f fa 



(82) 



Since the condition leading to this result is j^-jj 3> 1, the powers of (x^f) i n equations ([51]). 
will be smaller than one, and thus not an enhancing factor. From (|82p it appears that in 
the scenario of modified dispersion relations, the (nearly) folded configurations, singled out by the 
condition l + cos(#j) < Xfcf ^' are en hanced with respect to the other configurations, but not as 
much as in the modified vacuum scenario. In any case, the result shows a non-trivial Oj dependence 
(recall that we consider a realistic squeezing limit where k\ is small but non-zero). Furthermore, 
we see that the higher the number of derivatives, the more the suppression because of the factor 
(li") 1+mm ^ n ' m '^ : a § a i n ) the cubic couplings ([22]) and (pl3"|) that we have discussed as detailed examples 
turn out to be leading. 

thus, once more, the possibility 



Equation (|8ip also includes enhancement factors such as 



Ml 



HA 



or not of actual overall enhancements depends on the values of the scales and quantities at play. Let 
us again consider the same example as before and the couplings (|22p and (|43p . to give some pointers 



about the possibilities. Recall that we take e 



I/' 



10" 



H ~ 10" 5 M D 



compatible with the 



results of WMAP, and A ~ 10 M Planck (supersymmetric GUT scale). Then, if we can only probe 



values down to ^ > 10 2 , 



we are in the case > 1 and we find the bound, roughly, 



mv, mdr ^ 
I min cub r~^> 



10-^(1 -n s ) 
10(1 — n s ) mv 

10^+T (-77-) ~ mdr 



ifl+co S (^)>fi 
ifl+co S (^)<fi 



(83) 



We see that for |^ ~ 10 -2 the folded configurations, both in the modified vacuum and modified 
dispersion relations scenarios, could be ten time enhanced with respect to the standard result, whereas 
the other configurations would present only a moderate enhancement compared to the result in the 
standard scenario. 
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In the case of the higher-derivative interaction (|43p . for a coupling g ~ 10 2 , we find, see (j-52D . 
54")) and the comment at the end of section 14.3. 2\ after (|72p , 



mv, mdr ^ 
I quart dcr 



10-^(1 -n s ) 
fci 



10 J 



10 d 



-l , 



- 1 (-^)' c + 1 x (l-n s ) mdr 



ifl+cos(^)>fi 



ifl+cos(^)<f^ 



For |±- ~ 1(T 2 , 



(84) 



and 



are of the same magnitude, which again confirms that in the mod- 



ified scenarios the different a -1 scaling of the higher-derivative interactions can compensate for the 
stronger suppression by higher powers of ^ and make the squeezing limit sensitive to higher-derivative 
interactions. 



5.2 Dependence on k± 

We discuss now the dependence of the corrections to the bispectrum in the squeezing limit on the 
squeezed wavenumber k\. Recall that although fik s depends on the scale k$ in a way that we cannot 
fully determine with our phenomenological approach, it does not depend on k±, thus the /ci-dependence 
of the leading correction to the bispectrum is not affected by our ignorance of the scale dependence 
of (3 ks . 

In the following, we will rewrite equations (|6"6"|) - ([6"T|) and (JTT)) - (fT2"j) highlighting the dependence on 
k\. Again, we find that the field theoretical results for the modified scenarios have two different kinds 
of behaviour, depending on the magnitude of the scale of new physics as compared to the sensitivity 
of the observations. We will show that the results are different from the one in the standard scenario, 
which grows as ~ fcj -3 in the squeezing limit [IHBj, and from the one obtained with the folded template 
proposed in [9] to model modified vacuum scenarios, which yields a bispectrum growing as ~ k^ 2 [HE]- 

Cases <C 1 and f 1 -^ < 1 

In this case, equations (|66l)-(l67D and <^-$T2§ yield 

(85) 

with D = \k s r] c \ or respectively for modified vacuum or modified dispersion relations. 

2 — |— YTl 1 II ^ Tb 7Tb S ^ ~\~ Tt " 

The correction scales as k 1 t , = {0, 1}. We thus find that higher-derivative 

couplings (larger values for n, m, s) grow more slowly. The fastest growth for small k\ is of the form 

A M 2 / H D\ n+m+s 
S^MkM^^^sl^^^^) PAh)Pjk s ), (86) 

that is, a local shape fc^~ 3 , and can be found for couplings such as the minimal coupling cubic, with 
min(n, m, s) = — 1 (recall that this value occurs due to the combination (d)~ 2 ip, therefore it is always 
accompanied by nf = 1). 

Also in these cases, however, although it has a local form, the result is different from the one in 
the standard scenario because non-Gaussianities can be enhanced. 
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Cases \kiTj c \ > 1 and » 1 

In this case, equations (|66|) - ([67|) and ([7T|) -([72 |) yield, at leading order, 



' ks \ — 1— mm(n,m,s) 



Ml 



Hks_ 
A ki 



V2e HA 

n+m+s 



\ A hi) 



n + m. + s 

'Hks\ K+1 



mdr 



if 1 + C0S&-) > f £ 



ifl+cos(^)<f ^ 



(87) 



Except for the case of folded configurations in modified vacuum scenarios, to which the same consid- 
erations regarding equation (|85|) apply, it appears that higher-derivative couplings (larger values for 
n, m, s) would lead to a more pronounced growth for small k\ as its exponent will be more negative. 
However, all these contributions are very suppressed because of the factors (# ) n + m+s ; (# ) «+i _ ^he 
most interesting cases are thus those that are least suppressed, like the minimal cubic coupling that 
we studied as a detailed example, where n + m + s = 1 and min(n, m,s) = —1, see (|60|) . in which case 



M 

-I/O I Planck 

£ |Pfcsl — ^ — 



r # m 

A fcj 



ifl+cos(^)>fi 



Af) K+1 m dr 



if l + cos(^) < 



H_J_ 

A xx 



\P st {k x )P st {k s ). 



We see that in this case the form of the bispectrum in the squeezing limit is i) not local, with a scaling 
fcr/ 4 , and with amplitude not too severely suppressed for non- folded configurations, %%) enhanced (see 



also section I5.ip . with a non-local shape for (nearly) folded configurations in the case of modified 
dispersion relations, and Hi) with a local shape and greater enhancement for folded configurations in 
the modified vacuum case. Remarkably, the scaling with k\ of the contribution from (nearly) folded 
configurations in the case of modified dispersion relation captures the power k of the first momentum 
correction to the standard dispersion relation. 



5.3 Predictions from the different scenarios 

We have seen that, most interestingly, in modified scenarios the leading features of the bispectrum 
in the squeezing limit depend on very general aspects of the theory (exit from horizon, particle pro- 
duction, interference). These cannot determine the bispectrum in the squeezing limit in all details, 
but even adopting a wide-range phenomenological approach, we have been able to determine both 
the leading dependence on the squeezed scale k\ and the presence of and bounds on enhancements of 
non-Gaussianities. 

We have seen that these features depend on the values (less or more than 1) of the quantities 
l&i^cli Ysli- We would like now to ask ourselves which kind of conditions or predictions the different 
theoretical scenarios impose or make about the values of |fcir? c | an d f 3 ^) and, thus, what features of 
the bispectrum would be realized in the squeezing limit. 

• EBA modified vacuum. In the Effective Boundary Action approach, the boundary condition for the 
fields is fixed at a time rj c independently of the modes k ("beginning of inflation"). In that scenario, 
it is then assumed that all modes probed by observations were subhorizon at rj c , because in that 
case one can fix initial conditions with more confidence. Therefore, since also k± is an observable 
mode, it must be \k\r} c \ > 1 (subhorizon mode at r? c ). 
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Therefore, if the EBA is the correct way to model the effects of very high energy physics from a 
low-energy point of view encoding them in an initial state for the perturbations, then the leading 
correction to the bispectrum for non-folded configurations would grow at large scales as fcj~ 4 , but with 
only moderately enhanced non-Gaussianities, whereas the contribution from folded configurations 
would be more enhanced, but growing only as k^ 3 , see ([80]). ([83]) . ([87]) . (|88|) . 

• NPHS modified vacuum. In the NPHS approach, instead, the boundary conditions is imposed at 
the time given by \kr] c (k)\ = jj 3> 1. In this case, the picture is that different modes are created 
out of the vacuum at different times in a certain condition that fixes the vacuum of the effective 
theory at scales lower than A. When computing the three-point function, therefore, we cannot trace 
the interaction of the modes backward in time for a time earlier that the "creation time" for the 
largest of the considered modes: the prior evolution (and our ignorance about it) is encoded in the 
initial state. Thus, in the case of the squeezed limit, the earliest time possible is related to the large 
momentum scale & max = m&x(k l 23 ) ~ k$, and not to k\, so that f^-jy can be both larger or smaller 
than 1, depending on the values of the scales at play in the specific model. 

In particular, the leading correction to the bispectrum would grow at large scales as k^ 4 for non- 
folded configurations, only provided that the scale of new physics is such that A > ^H, where jj^ is 
the smallest ratio of scales that can be probed, see ()87|) . (|88p . For folded configurations, and when 
A < ^H, the contribution would instead be of the local type, growing at most as k^ 3 . 

An interesting outcome of this result, is that in case we could detect one or the other behaviour, say 
in the halo bias, and separate it from other sources of non-Gaussianities, we could obtain information 
on the magnitude of A, given our knowledge of the observation sensitivity 

• Modified dispersion relations. In the scenarios with modified dispersion relations violating WKB at 
early times, in case of a Lorentz-breaking scale A such that > 1, where ^ is the smallest ratio 

that can be probed, there would be contributions to the bispectrum growing as fast as ki 3 ~ N with N 
only bound by the specific model, and in principle also possibly large, due to the presence of higher- 
derivative Lorentz-breaking couplings, see (|87[) . However, we have found that those contributions 
are severely suppressed by powers (jr) , M > N. The least suppressed contributions (which could 
actually be enhanced) show at leading order a growth as k^ A for non-unfolded configurations, and 
slower powers k^ r , with r fractional, for the (nearly) folded ones, see (|87|) . (|88|) . Remarkably, in the 
latter case, the exponent of k\ captures, being directly related to it, the momentum power of the 
lowest correction to the standard dispersion relation. 

On the other hand, for A < jr^H, the leading contribution would be of the local type growing at 
most as k^ 3 , see ([86]) . Once again, if we could observe one or the other behaviour, we would obtain 
information on the magnitude of the Lorentz-breaking scale from our knowledge of the probed k\ . 



6 Comments on the halo bias 

We will here briefly comment on the implications of our results for the halo bias, leaving a more 
detailed analysis for future research. The halo bias depends on the bispectrum and the spectrum 
as [7] ([8] has a similar but not identical formula) 



Ab h (k u R) 5 C 1 



D(z)M R (k 1 ) 8vr 2 (j 2 



R JO 



dk 2 klM R {k 2 ) 



W fit + U> + 2 fab{ ) ^ ^ ±|± 2hh ^ . (89) 
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where a\ is the variance of the dark matter density perturbations smoothed on a scale of Lagrangian 
radius R, 5 C is the critical threshold for the collapse of a spherical object (for a matter-dominated 
universe 5 C = 1.686), D(z) is the linear growth factor normalized to be D(z) = (1 + z)~ l during 
matter domination, P and B are the two-point correlatoi^l and bispectrum for £, and Mr is the linear 
relation between the dark matter density perturbations smoothed on a scale R and the primordial 
curvature perturbation: 

2k 2 

M R (k) = ^Jj2 T (k)WR(k) , (90) 

where Q m is the present time fractional density of matter. Finally, Hq is the present Hubble rate, T is 
the transfer function normalized to one on large scales and Wr is the filter function with characteristic 
scale R. 

The halo bias depends on both the primordial bispectrum and the subsequent gravitational pro- 
cessing which generates additional non-Gaussianities. These General Relativity corrections have been 
studied for example in [7] and it has been shown that they give a typical contribution to the halo bias, 
which could permit to distinguish them. We will therefore comment only the primordial contributions. 

The behaviour of the halo bias at large scales has been computed for the local, equilateral and 
folded templates for the primordial non-Gaussianities in [7] (see also [8]), finding that the bias goes 
approximately as k± 2 with the local template, as a constant in k\ with the equilateral, and as k^ 1 
with the folded (the folded template goes indeed as k^ 2 ). 

The latter result is the one that can be compared with ours in the case of the modified vacuum 
approach, since the template was proposed as an approximation to model the bispectrum for modified 
vacuum scenarios [9]. In fact, we have shown that this leads to the wrong result in the case of the 
squeezing limit: the template does not depend on the additional scale(s) i] c ,A and thus its squeezing 
limit is different than the one that actually occurs in the rigorous field theory computation. We 
have indeed shown in sections 14.1.21 14.2.11 and 14.3.21 that for realistic limits (in view of realistic 
observations), the leading contribution to the bispectrum grows as k^ 4 or k^ 3 , depending on the 
compared magnitudes of the scale of new physics and the sensitivity of the observation (smallest 
squeezed momentum scale), so that the halo bias would approximately grow either as k^ ?J or k± 2 , and 
not as kr[ . 

Similarly, in the scenario of modified dispersion relations violating WKB at early times, the be- 
haviour of the bispectrum in the squeezed limit, see sections 14.1.31 14.2.21 14.3.2) would also lead to a 
bias that approximately grows either as A;^ 3 or k^ 2 for small k\ for non-folded configurations, again 
depending on the relative magnitudes of the high-energy scales and the sensitivity of the observa- 
tions. (Nearly) folded configurations would instead show different scalings with k\, determined by the 
exponent of the first correction to the standard dispersion relation for small momenta. 

In general, a feature relevant for the observations of the bias is that the high energy modifications 
could enhance the single- field non-Gaussianities in the squeezing limit, compared to the practically 
undetectable level that is present in the standard scenarios of single field models with Bunch-Davies 
vacuum and standard kinetic terms. However, the final magnitude of the non-Gaussianities depends 
as well on the specific model, which determines the precise form and magnitude of the coefficients fik s 
entering the equations. 

7 Final discussion and conclusion 

In this article we have analized the squeezing limit of the bispectrum in single-fields models of infla- 
tion when modifications of the theory at very high energy are present. Due to the large number of 

20 In the literature on the halo bias, the two-point correlator is often called the spectrum [7][8], not to be confused with 
the quantity V = £^P, also called spectrum, usual in the CMBR literature. 
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inflationary models available, it is useful to look at general features that can be obtained when these 
modifications are rather generic and simple. In particular, we have considered 1) the presence of terms 
with higher derivatives, which are not suppressed any more beyond a certain momentum scale and 
modify the dispersion relations of the fields, and 2) the possibility of initial condition for the solution 
of the field equation that parametrize in a simple and generic way our ignorance on the physics at 
very high energies. The leading contribution to the squeezing limit appears dominated by very general 
features of the scenarios (particle creation at early times, interference) and thus it is possible to obtain 
general results and bounds. 

Those that we have found are interesting in two respect: on the one hand, they show that non- 
Gaussianities could be enhanced by high energy modifications of the theory, which would improve the 
possibility of detection of those phenomena. On the other, the behaviour of the bispectrum at large 
scale in these scenarios is distinctive and differs from the one that has been found in the standard 
scenarios, with Bunch-Davies vacuum and standard dispersion relation, and also different from that 
found in the previous literature using the approximate template proposed in [9] for the modified 
vacuum approach. Particle creation at early times and interference effects explain the differences 
with the standard result. Instead, the differences with the result obtained using the template for 
the modified vacuum approach arise principally because the approximate template does not actually 
depend on the scale of the new physics, while the bispectrum calculated from the field theory does, 
and therefore taking the squeezing limit in the template cannot account for the interplay between the 
scale of the small mode k\ and that of the new physics when taking the limit itself. 

Our results show that it could be possible to obtain various different pieces of information about 
the features of the high energy physics theory from large scale observations such as those of the halo 
bias, if the signatures we have found were detected. In fact, the dependence of the bispectrum on the 
small wavenumber k\ in the realistically probed squeezing limit changes in relation with the magnitude 
of the new physics scale: if able to detect certain dependences of the bispectrum, one can place bounds 
on the new scale of physics. 

Of course, other sources could give origin to sizable non-Gaussianities of non-local shape, for 
example in multi-field models, and thus an interesting question for future research would then be if 
the effects due to high-energy modifications of the theory that we have considered can be distinguished 
from those. For example, the dependence on the small mode k\ of the halo bias could be different, 
or the enhancement factors could be absent or of different magnitude. We leave this as an interesting 
outlook for future research. 
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A Appendices 

A.l Solution of the field equation in the case of modified dispersion relations 

We parametrize the modifications to the dispersion relation via a generic function F {^krj): 

u,(rj, k) = a(r,)u, pbya (p) = a{r,)p F (-£) = k F (^-kr^j , F(x ->• 0) -)■ 1 , ^ « 1 , (91) 

with the only request that the WKB conditions be violated for a short period of time at early times. 
The generic shape of such a dispersion relation is in figure [TJ The approximate solution of the equation 
of motion ([5]) in these cases is given in equation ()10p . where, see [16] : 
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in regions I and III, in terms of the convenient variable ^rkr] = yk- 
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V 2k u(y k ) 
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(92) 
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(94) 



in region II, the partial solution U\ )2 {r],k) can be found either by solving in details the specific 
equation, if one has a favourite model, or, more generally, by expanding the frequency around 
the minimum wo at 77 = rj m as 



1 



W 



1 1 

aw, ~ 9rr9 

phys ^jZJ-[2 



(95) 



so that 



fk(v)\ 



y^BiUi(ri,k) = B\ W(iKrj m ,a,2iKr]) + B 2 W(-iKr] m ,a,-2inrj) 
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^Jm 2 - 4oj 2 - 4H 2 K 2 r)l 



H 



2H 



(96) 



where W(a,b,z) is the Whittacker function. As explained in [16], this approximation is well- 
justified also because backreaction constraints the interval [?7i,?7n] to be very small (A < 1), so 
that rj ~ r\ m for 77, rj m € iy n ]. In any case, it actually turns out that the details of the solution 
in region II are not important for what concerns the leading contribution to the bispectrum |16j . 



By asking for the continuity of the function and its first derivative, and imposing the Wronskian 
condition W{/, /*} = — i to have the standard commutation relations in the quantum theory, we 
obtain in full generality 



oik 



W{Ui,U 2 } 
W{B 1 Ui+B 2 U 2 ,u 2 } 
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(97) 
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(99) 



>ln 



where W is the Wronskian. We choose ^ = 1 picking up the usual adiabatic vacuum. The Wronskian 
condition also imposes a 2 — fif. = 1. By expanding for small A around rji, we obtain equation ([12]) in 
the text. 

Finally, the parameter signalling the WKB violation is |16| 



1 



— + 3— 

4w 2 8w 3 urj 2 



+ 2w 4w 2 2U + AU 2 ' 



(100) 
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A. 2 Formulas for Whightman functions and their time derivatives 

The bispectrum is computed at a time rj when all momenta have exit the horizon, therefore the 
Whightman functions have the general form 

r?^o ^ a(r/) a(rf) 

We list here the explicit formulas obtained substituting in the above equations the field modes 
functions fk for Q in the various scenarios we consider. Due to the presence of positive- and negative- 
energy solutions, the Whightman function will also have a positive- and negative-energy, which we 
indicate with G^. 

• Standard scenario and modified vacuum. Prom equations ([7]), (|14p 



G±(0 )I) ') = ±^(lT^' d^^n') = ±S^^k 2 v >e ± ^', (102) 



• Modified dispersion relations. From equation (|10p , we find that in region IV the Whightman function 
is the same as the one in equation (|102p . while in region III it is obtained from (|92p . Thus, 

(p uyri ) \{ Tt'J 111 re S lon IV 

The derivative of the Whightman function can be written in a form valid both for regions IV and 
III, given the limiting behaviour of and U at early times, so that 

<p z a {V ) [i m region IV 

(104) 

where the suffix M means that the positive-energy solution is associated with 7*,x*> while the 
negative-energy with 7, x- Finally, T~L is the conformal Hubble scale. 

A. 3 Time integral in bispectrum computations with modified dispersion relations 

When computing the bispectrum, we have encountered various Laplace integrals of the form 

M w) = [" ' dy'y'^e i ^ {y ' ) h{x 1 ^9 j: y'), (105) 



see equations (|35|) . ([53]) . (f70|) . where w and h(xi,xs,0j,y') change in the different integrals and the 
variables are defined in ([34"|) . 

In particular, in the case of the minimal coupling interaction (|22p . the integral present in the 
correction of order j3]~ to the bispectrum, see equation ([55]) . has the form of (|105p with w = 1 and 
h({x},9 jt i/) 9{{x},e 3 ,y'), where [E] 



g({x h ^}, Xj ,y',m) = ]J j*(x h , y')i{x 3 ,y') e*f (^ s M+ {-i) m s^)) (1Q6) 
with 7 defined in (|104p and 

4u;{ Xi yy 8w(x i y)' :i u){Xiy)y 2 \j(x,y') if m = 1 
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\y\ < i 


y = yn 
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Tabic 1: Behaviour of functions entering the integrands in equations ([55)1 , (|53p . ([70]) close to the boundary 
points (limits of integration) j/ n ~ -1 and j/ ~ -C -1, and close to the possible stationary point (s) 
y* in the squeezing limit x\ -C 1, xs — 1- 



In the squeezing limit, for m = 1 and j = 2, 3, looking at (|93|) and (|104|) . we obtain 

g{y) = #(xi<l,x. „ 3 ~x s ~l,y,l) ~ -i |7(z s ,y)| 2 e** °( X1 'T ) ~ -i| 7 (l,y)| 2 (108) 



In the case, instead, of the higher derivative interaction in equation (j43j) . the integral appearing 
in equation ([53]) has tu = 3 and h({x},9j,y') — > q({x},8j,y'), where [T6] : 

J* X i + l X i+2 X i+ 1 Xj 

+ ^ ± ^X*(^ + Mx ] h*(x j+1 ))e M* , (109) 

Xj / 

and 52 has been defined in equation (jl07j) and x i n ()103j) . We write explicitly only the dependence 
on x of 7, x to avoid cluttering the formula. The labels of the x's are defined modulo 3, here. 
In the squeezing limit, for j = 2, 3, looking at ([93]) and (|104p . 

q(9j,xi,y) = c7(xi<l,x j=2 3 ~x s ^l,y) ^ -6i\j*(l,y)\ 2 +2i — - - 4ix*(xi, y) cos(6j) (110) 

u [.*■} y ) 

For ^xi 3> 1, the integrals like (|105|) are dominated by the contribution of boundary and stationary 
points, where the oscillations due to the large phase of the integrand are reduced. It is therefore 
important to study the behavior at leading order of the integrands in ([35]) . ([53]) . ([70]) close to the 
boundary points (limits of integration) y n ~ — 1 and y ~ <C — 1 and to the possible stationary 
point(s) {y*}. In table [1] we report the leading and next-to-leading behaviour. 

In the table, we include also the function F(x, y'), defined in (|9ip . and now written in terms of the 
variables (|34p . From its definition (|9ip . it appears that the expansion in powers of ^ = — y for \y\ <C 1 
is in fact the expansion for the effective frequency from the effective action. In particular, k+1 is thus 
the power of the first correction to the standard physical dispersion relation (model dependent): 

u ph Jp) ~ p(l + (f) K FW + ...). (Ill) 
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In expanding g(y), q(8j,x\,y) - see (|108|) . (jllOp - we have taken into account the WKB conditions. 
The powers z/*, //*,r* are model dependent. We wish to be generic, so k, v*, t* can be fractional or 
integer^!. Then, since F = Uv ^ a is a real function, its expansion must be clearly in powers of \y\, \y—y*\, 
with real coefficients, for generic v*, k. Similarly, it must be so for v e , (x, y)A = — (p + w p h ys (p)), which 
is a real function as well (in table [1] we have been very generic concerning its expansion around y*). 
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